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Precision measurements are performed on highly excited vibrational quantum states of molecular 
hydrogen. The v = 12, J — 0 — 3 rovibrational levels of H 2 lying only 2000 cm“^ below 

the first dissociation limit, were populated by photodissociation of H 2 S and their level energies were 
accurately determined by two-photon Doppler-free spectroscopy. A comparison between the exper¬ 
imental results on u = 12 level energies with the best ab initio calculations shows good agreement, 
where the present experimental accuracy of 3.5 x 10“^ cm“^ is more precise than theory, hence 
providing a gateway to further test theoretical advances in this benchmark quantum system. 
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Quantum chemistry started as a pioneering applica¬ 
tion of quantum mechanics, when Heitler and London 
explained in 1927 the existence of an attractive bond¬ 
ing state between the two hydrogen atoms to form the 
ground state of molecular hydrogen.® Since then, H 2 has 
been a canonical test system, where innovations in ab 
initio theory, most notably by Kolos and WolniewicJ^ 
and experimental measurement techniques, most notably 
by Herzberg and coworker^^l® have mutually stimulated 
each other as a driving force toward the highest ac¬ 
curacies.® Similar highly-accurate comparisons between 
theory and experiments are also actively pursued in 
molecular hydrogen ions and HD+, where accura¬ 
cies achieved in calculations are better, owing to their 
simpler three-body configuratioiPMi^. This simplicity on 
the other hand, makes the one-electron ion system to 
be atypical of molecules as it does not feature electron 
correlations, which form an essential ingredient in the 
quantum chemistry of molecules and condensed matter 
systems for which treatment of H 2 serves as a benchmark. 

In recent years, great progress has been made in the 
quantum chemical calculations of the energy level struc¬ 
ture of the H 2 ground electronic state. For exam¬ 

ple, the most recent calculation of the chemical bonding 
energy or dissociation limit Dq of the ground electronic 
state has an accuracy at the 10“^ cm“^ level, while nearly 
equally accurate binding energies were calcul ated for the 
entire manifold of rovibrational states.^^^^ To achieve 
these accuracies, the Born-Oppenheimer (BO) potential 
energy data points were calculated to 10“® cm“^ ac- 
curac resulting in BO-level energies of 10 ® cm ^ 
accuracy. Calculations of adiabatic corrections to the 
BO energjISl yield accuracies of 10“® cm“^, while non- 
adiabatic energy corrections were determined at uncer¬ 
tainties at the ~ 10“"* cm“^ level.^ At such level of pre¬ 
cision the relativistic and radiative or quantum electro¬ 
dynamic (QED) effects, amounting to several 0.1 cm“^, 
need to be accurately accounted for, as was accomplished 
in the novel approach by Komasa et 

Direct purely vibrational transitions in the ground 
electronic state of H 2 are extremely weak due to their 


quadrupole nature. They were predictecP^ and subse¬ 
quently observed by Herzberg.E3 The studies of Bragg 
et using a Fourier-Transform spectrometer com¬ 

bined with high-pressure absorption cells exhibiting ef¬ 
fective path lengths of 400 m, enabled the measurement 
of the fundamental and overtone vibrational splittings up 
to the (4-0) band. This was further extended to the (5-0) 
band in laser-based studieJ^^ using multipass cells with 
effective absorption paths of 20 m. Recent laser-based 
measurements have improved the experimental accura¬ 
cies for the fund amen tal (1-0) vibrational splitting using 
molecular beam^^nHH reaching uncertainties of 1 x 10“^ 
cm“^. Sensitive cavity-ring down laser spectroscopy was 
applied to the (2-0) overtone bancP^ at uncertainties of 
1 X 10“^ cm“^, as well as on the (3-0) overton e ban d 
achieving ^ 10“® — 10““^ cm“^ uncertainty levels 

It is unlikely that the direct quadrupole excitations 
can be extended up to the highest vibrational quanta, 
since for example the (12-0) transition in H 2 is some six 
orders of magnitudes weaker than the fundamental (1-0) 
band.l^ Thus far, the H 2 level energies obtained fro m the 
calculational framework of Pachucki and coworker j^^O^ I 
are in excellent agreement with the experimental de¬ 
termination of the dissociation energy,^ the fundamen¬ 
tal vibrational splitting,!^ the (2-0) and (3-0) overtone 
band transitions, I22H211 as well the level energies of the 
rotational series of the lowest vibrational state.^^ These 
highly-accurate measurements test the quantum chemi¬ 
cal calculations of the low vibrational quantum numbers 
v = t)-3. 

The present experiment seeks to test the calculations 
of H 2 binding energies for vibrational levels v = 6 — 12, 
where theoretical calculationJi^ exhibit the largest un¬ 
certainties. We follow up on the experimental findings of 
Steadman and Baer,l^ who first produced highly vibra- 
tionally excited hydrogen (II 2 *) from the photolysis of 
H 2 S. Besides the main photo-dissociation product chan¬ 
nel yielding SH molecules, there exists an energetically al¬ 
lowed channel producing H 2 * under non-equilibrium con¬ 
ditions, upon absorption of two UV photons in the H 2 S 
molecule. They performed a single-color excitation ex- 
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FIG. 1: Schematic of the measurement setup showing the 
various lasers used in the photolysis of H 2 S (292 nm), for 
performing the EF — X spectroscopy step in vibrationally- 
hot H 2 (291 - 293 nm), and for inducing the ionization step 
(202 nm). The probe laser is split and interferometrically 
aligned in counter-propagating fashion to obtain Doppler-free 
two-photon excitation. Resonantly produced H"'"-ions are de¬ 
tected after a time-of-flight (TOF) region for resolving the ion 
masses. 


periment, in which the UV-pulses caused two-photon dis¬ 
sociation, two-photon excitation in the H 2 EF—X(y\ v") 
band, and multi-photon ionization/dissociation produc¬ 
ing ions. These sequential processes may involve 
some 6 or 7 UV photons near 290 nm depending on the 
excitation pathway. 

In our present study, for which a setup is displayed 
in Fig. [2 we disentangle the laser excitation processes 
using up to three different tunable UV-laser systems to 
perform a controlled study aiming at spectroscopic accu¬ 
racy. In all cases, a powerful UV-laser pulse, obtained 
from a frequency doubled pulsed-dye laser and set at 
292 nm, is used for optimal photolysis of H 2 S and effi¬ 
cient production of H 2 *. The other laser pulses are time- 
delayed in steps of ^ 10 ns with respect to the photolysis 
laser to avoid temporal overlap, and therewith circum¬ 
vent broadening of the lines. Signal is obtained by moni¬ 
toring hydrogen ions, produced in the interaction region, 
accelerated by ion lenses and propagating over a time-of- 
flight (TOF) mass separation region onto a multichannel 
plate detection system, where the mass-resolved signal 
is recorded. In virtually all cases the signal is most 
prominent, although on some lines also signal is de¬ 
tected. 

As a first step of our investigation the output of a sec¬ 
ond frequency-doubled dye laser was spatially overlapped 
in the interaction region with the photolysis laser. A 
scan over the two-photon excitation region 68 200—68 800 
cm“^ region is displayed in Fig. The lines in this low- 
resolution spectrum exhibit broadening due to the laser 



FIG. 2: a) An overview spectrum of H 2 showing lines as¬ 
signed toF-A(l,10), F-A(2,ll), andF-A(3,12) bands. 
Detection is in the H'*’ channel. The most intense line, marked 
by an asterisk (*) and exhibiting off-scale intensity, is unas¬ 
signed. b) Potential energy curves of and elec¬ 

tronic states, including two-photon transitions in the F — X 
system probing high vibrational levels. 


linewidth, Doppler and ac-Stark effects, in the order of 
several cm“^, but are narrower than in the one-color 
study of Steadman and Baer.l^ At the high intensities 
used the ac-Stark effects also causes appreciable shifts to 
the transition frequencies, but the accuracy is sufficient 
for a tentative identification of the transitions. 

In the overview spectrum the two-photon lines in the 
(3,12), (2,11) and (1,10) bands of the F — X system were 
identified. Note that the numbering of vibrational levels 
in the FF double well system can be done separately for 
the inner F and outer F well or jointly for both wells; 
FF{5) then corresponds to F{3). These spectra and 
the much narrower lines measured in the Doppler-free 
spectra (see below) revealed many misassignments in the 
proposed assignments of Steadman and Baer ,1^ although 
the generic finding of H 2 * is confirmed. In particular 
the assignment of transitions to the F inner potential 
well could not be confirmed, and we conclude that it is 
doubtful that any F—X transitions were detected in their 
study.^ However, assignments of most transitions to the 
F outer well states are confirmed presently. Based on 
Franck-Condon arguments the vibrationally-excited lev¬ 
els in the ground state can be most easily excited 

to the outer well U-levels. The potential energy diagram 
of Fig. I^b) illustrates that the high vibrational ground 
state levels w = 10 — 12 have favorable Franck-Condon 
overlap with levels in the outer F well. In our recordings 
a number of lines remain unidentified, in particular the 
strongest resonance observed, which is marked with (*) 
in Fig.[^ These problems with assignments may be due 
to the fact that high(er) J-levels in the EF state, in par¬ 
ticular in the outer F well, are not known to sufficient 
precision to establish positive identifications. 

We choose the rotational series of Q transitions, i.e. 
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AJ = 0, in the F — A'(3,12) band as the main fo¬ 
cus for a precision study in H 2 with the aim of testing 
the recent quantum theoretical ab initio calculations!^ 
For the precision frequency measurements the beam of a 
UV laser, obtained from a narrowband injection-seeded 
traveling-wave pulsed dye amplifier (PDA), is split in two 
equal parts and interferometrically aligned in counter- 
propagating fashion, to perform two-photon Doppler-free 
spectroscopy. Most of the H+-signal is produced by the 
combination of the photolysis and spectroscopy lasers. 
In cases where the highest spectroscopic accuracy is de¬ 
sired, the power of the PDA was tuned down, and a third 
UV laser beam at 202 nm (frequency-tripled dye laser) 
is used to assist ionization from the ^"(3) state. To avoid 
broadening effects induced by this ionization laser, the 
pulse timing is again delayed by 10 ns with respect to 
the spectroscopy laser. In the expectation that possibly 
auto-ionizing resonances could help to increase signal the 
third laser beam was scanned over the window 202-206 
nm, but no such resonances were found. 

The two-photon Q-branch lines in the F — A(3,12) 
band were recorded under the improved high-resolution 
spectroscopic conditions using the narrowband PDA laser 
system with a bandwidth of < 100 MHz. The PDA 
system has been described recently, including the devel¬ 
opment and characterization of a novel frequency-chirp 
analysis procedure to identify possible offsets between the 
continuous-wave seed frequency and the pulsed output of 
the PDA. This allows for accuracies of 0.001 cm“^ to be 
achieved.!^ In Fig.j^a recording of the two-photon Q(l) 
transition in the F—X{3, 12) band is displayed, measured 
at different intensities of the spectroscopy laser. 

For absolute frequency calibration of the cw-seed laser, 
hyperfine-resolved saturation spectra of I 2 are simulta¬ 
neously recorded, using a split-off portion of the cw-seed 
radiation, and compared with the I 2 database.!^ At the 
same time, transmission peaks from a stabilized etalon 
with a free spectral range (FSR) of 150.01 (1) MHz pro¬ 
vide a relative frequency calibration, assisting to bridge 
the distance between H 2 and I 2 resonances.!^ The contri¬ 
bution to the error budget for the H 2 transition frequen¬ 
cies from the laser scan nonlinearity, etalon calibration, 
and the I 2 calibration are estimated to be as small as 2 
MHz. The chirp-induced frequency offset between cw- 
seed laser and the pulsed output of the PDA system is 
determined and corrected for in the final frequency, yield¬ 
ing typical values of 11 (3) MHz. Note that final correc¬ 
tions must be multiplied by a factor of 4, to account for 
harmonic up-conversion and two-photon excitation. 

The major source of line broadening and uncertainty 
of the transition frequencies is the ac-Stark effect. Since 
the transitions are weak and the population density of 
H 2 low, it was necessary to use relatively high probe 
laser intensities, causing substantial ac-Stark broaden¬ 
ing. With a laser bandwidth of ~ 100 MHz, the UV 
two-photon line-shapes are expected on the order of 
^ 200 MHz, but even at the lowest laser intensities 
of the combined counter-propagating beams the effec- 
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FIG. 3: (a) Measurements of the Q(l) two-photon Doppler- 
free transition in the F — A(3,12) band, (b) Plot of the 
peak positions at different intensities; the intercept at zero 
intensity (evaluated for both linear and quadratic fits), yields 
the (ac-Stark-free) transition frequency. 

five line widths were ~ 1 GHz. Fig. displays record¬ 
ings of the Q(l) two-photon transition at different in¬ 
tensities in the interaction zone. At the lower inten¬ 
sities used the resonance frequency is determined from 
a Gaussian fit to the recorded spectrum, while at the 
highest intensities the resulting asymmetric lines shapes 
were fitted to skewed Gaussian functions of the form 
fsG = (expai^/2)/(-\/^) • [1 -t- erf(a;/-\/2)]. The result¬ 
ing peak positions are plotted against the relative UV 
intensity as shown in Fig. [^b). To account for the in¬ 
creased widths and skew, the peak position uncertainty 
in the high intensity scans are larger. In view of the long 
extrapolation to zero intensity fits were made for a lin¬ 
ear and a quadratic fit. Results from these approaches 
were averaged and the uncertainties were conservatively 
estimated to cover the values when differing. In Table |l] 
the thus determined two-photon transition frequencies 
for the Q(0)-Q(3) lines in the F — X{3, 12) band are listed 
with accuracies of typically 100 MHz or 0.0035 cm“^ de¬ 
riving fully from this ac-Stark analysis. 

In order to make a comparison with the quantum 
chemical ab initio calculations these values for the ac¬ 
curate transition frequencies must be translated into val¬ 
ues of level energies in H 2 X^E+(?; = 12, J). This can 
be accomplished by subtracting transition energies from 
the accurate level energies of the upper F"(3) states from 
Bailly et al.,!^ which were determined via combined laser- 
based and Fourier-transform emission experiments yield¬ 
ing values accurate to within 0.001 cm“^ for J < 3. So 
this transformation does not add to the uncertainty of 
H 2 * level energies, which are listed in Table [h] as ex- 
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TABLE I: Measured two-photon transition frequencies in the 
F — X{3, 12) band, equivalent to the EF — X(5,12) band. 

Line Frequency 
Q(0) 68 476.0459(35) 

Q(l) 68 446.2834(35) 

Q(2) 68 387.6623(35) 

Q(3) 68 302.0250(35) 

S(l) 68 505.518(7) 


TABLE II: Level energies (Eexp) of ro-vibrational states 
(u = 12, J = 0 — 3), experimentally determined from 
the measured F — X{3, 12) two-photon Q-transitions and the 
accurate F{3) levels from Bailly et alf^ Also indicated are 
the theory values (Lithe) of Komasa et aZ.,^ and the difference 
AEexp-the between the experimental and theoretical values. 
The uncertainties in the differences represent the combined 
uncertainties from experiment and theory taken in quadra¬ 
ture. 


J" 

K 



-E'the 


^-^exp —the 

0 

34 302.1823 

(35) 

34 

302.1741 

(47) 

0.008 

(6) 

1 

34 343.8531 

(35) 

34 

343.848 3 

(46) 

0.005 

(6) 

2 

34426.2216 

(35) 

34 

426.217 9 

(46) 

0.004 

(6) 

3 

34 547.3362 

(35) 

34 

547.333 2 

(45) 

0.003 

(6) 


perimental values (i?exp) for v = 12. This procedure of 
combining differences between F(3) level energie^^ and 
X{12) level energieJ^^ provides a final and unambiguous 
assignment of the observed lines. The measurement of a 
weaker transition S(l) at in the same F — X{3, 12) band 
(A J = 2 for S transitions) allows for an additional consis¬ 
tency check of the line identification. The combination 
difference between S(l) and Q(3) transitions delivers a 
level splitting between J = 1 — 3 in X(12) amounting to 
203.494 cm“^, in agreement, at 0.008 (8) cm“^, from the 
ab initio values.l^In addition, the combination difference 
between S(l) and Q(l) lines is in reasonable agreement 
with the splitting between J = 1 — 3 in F(3) within 
0.010(8) cm-h 

In Table [ll] a final comparison is made between exper¬ 
imentally determined level energies (i?exp) and results 
from the quantum chemical calculations (Lithe) by Ko¬ 
masa et aZ.E^ For three of the four levels determined the 
difference between the experimental and theoretical val¬ 
ues, ALlexp-thej is Smaller than their combined uncer¬ 
tainties, while the J" = 0 level deviated by 1.3cr. This 
allows for the main conclusion of the present study that 
good agreement is established between experiment and 
the quantum chemical calculations for v = 12 levels in 
the electronic ground state in H 2 . This is the highest 
vibrational level subjected to a precision test below 0.01 
cm“^ for the H 2 molecule. 

In Fig. 1^ the contributing energy corrections to the ac¬ 
curate BO-energieJi^ are plotted as a function of v, as 



FIG. 4: Plot of corrections to the ah initio Born-Oppenheimer 
energies of the H 2 X^E^ ground electronic state^as a func¬ 
tion of vibrational quantum number v (in all cases J = 0). 
Adiabatic (ad), nonadiabatic (nad), and relativistic and ra¬ 
diative (rel+QED) corrections are plotted with reference to 
the scale on the left. The total uncertainty of the calculations 
(unc) is plotted with reference to the scale on the right. 


well as the final uncertainties in the theoretical binding 
energies.^ It shows that the calculations are least accu¬ 
rate for the range n = 6 — 12. The present test on u = 12 
covers this weakest spot in the calculations. Specifically 
the non-adiabatic contributions to the relativistic and 
QED energies (the so-called recoil corrections), as well as 
higher order non-adiabatic, a^-QED contributions, and 
the leading term in a^-QED, are mainly held responsi¬ 
ble for possible deviations with theory.E^ The presently 
determined experimental binding energies, being slightly 
more accurate than theory, may serve to stimulate cal¬ 
culations of such terms. 

As an outlook from the experimental perspective, ef¬ 
forts are underway to determine level energies for states 
up to the highest possible (y = 14) vibrational quantum 
states in H 2 . At least an order of magnitude improve¬ 
ment is foreseen in the experimental determination of 
the ?; = 10 — 14, J level energies, which should be achiev¬ 
able with an improved detection sensitivity, whereby the 
limiting factor of uncertainty, the ac-Stark broadening, 
can be suppressed. 

Note added during review. After submission of this 
manuscript a paper reporting recalculation of the non- 
adiabatic interaction was publishecPi^ improving their ac¬ 
curacy to the 10“^ cm“^ level. It was not specified in 
how far the newly calculated non-adiabatic wave func¬ 
tions affect the uncertainties of the QED and relativistic 
corrections. 
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tion for Fundamental Research of Matter (FOM) and by 
the Dutch Astrochemistry Network (NWO). We thank F. 
Ramirez (University of San Carlos, Cebu, Philippines) for 
assistance in the initial investigations. 
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